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We consider a hybrid single-electron transistor (SET) constituted by a gate-controlled normal-metal island 
(N) connected to two voltage-biased superconducting leads (S) by means of two tunnel junctions (S-I-N-l-S), 
operated as a turnstile. We show that the exchange of photons between this system and the high-temperature 
electromagnetic environment where it is embedded enhances Andreev reflection, thereby limiting the single¬ 
electron tunneling accuracy. 
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I. INTRODUCTION 


The experimental realization of a quantum electric current 
standard is one of the scientific and technological challenges 
of the present time. This is a key goal in metrology because 
it would lead to a modem definition of Ampere as well as to 
the most accurate comparison of the fundamental constants 
Rk = h/e 2 and Kj = 2e/hX Among the devices proposed un¬ 
til now (St; 7 the hybrid SINIS single-electron transistor (SET) 
depicted in Fig. Q] is one of the most interesting candidates^ 
Such a device is formed of a normal-metal (N) island joined 
to two superconducting (S) electrodes via two tunnel junc¬ 
tions with capacitances Cj for the source (S) and Cd for the 
drain (D). The entire structure is biased with a constant volt¬ 
age Vd — Vs = V. The amount of electric charge localized on 
the island is controlled using a gate potential V g (t), capaci- 
tively coupled to N by means of a gate with capacitance C g . 
Typically, the charging energy of the island Eq = e 2 /2CY, with 
Cy_ = C$ + Cd + C g , governs the tunneling processes in the 
SET, i.e., the system works in the Coulomb blockade regime. 
Additional energy filtering is provided by the two outer super¬ 
conductors which protect the device against unwanted tunnel¬ 
ing events. In this context, if the single-electron tunneling is 
the dominant process, a periodic V g (t ) signal with frequency / 
generates an electric current I through the SET which is equal 
to ef. In other words, the SET is a frequency-to-current con¬ 
verter. However, high-order tunneling events occur in addition 
to the single-particle ones. They limit the conversion accuracy 
of this electronic turnstile thereby acting as error sources. The 
main contribution to the total error is usually provided by elas¬ 
tic and inelastic cotunneling-^ as well as Andreev reflection 
and Cooper-pair cotunnelingJi- 12 From the theoretical point 
of view, it has been shown that all these processes can be elim¬ 
inated efficiently thereby reaching the metrological require¬ 
ments Nevertheless, in real experiments the achievement of 
the accuracy needed for the completion of the so-called quan¬ 
tum metrological triangle remains a difficult task. In partic¬ 
ular, a noise-induced residual Andreev tunneling current af¬ 
fects the I-V characteristic of the SET turnstile although the 
increase of the charging energy Ec, with respect to the gap 
parameter A of the superconductors, leads to a decrease of 
Andreev reflection probability^ Such a two-electron current 


may be due to the effect of the high-temperature electromag¬ 
netic environment the SINIS device is coupled with. The en¬ 
ergy provided by such an external thermal bath to the SET via 
the exchange of photons can promote tunneling of particles 
through the single junction^ - — In this paper, we show that, 
indeed, the environment-assisted Andreev reflection limits the 
turnstile accuracy, unless it is properly taken care of. 

source island drain 
Os Cd 


V d = V/2 



FIG. f. Hybrid S-I-N-I-S single-electron transistor (SET). The black 
parts stand for the insulating barriers of the tunnel junctions. 


II. ELECTRONIC TRANSPORT IN A SINIS TURNSTILE 

In the Coulomb blockade regime, the electronic transport 
in the SINIS device of Fig. [T] is determined by the charging 
energy Ec . For a symmetric device, C$ = Cjj = C, assuming 
that initially the excess electric charge localized on the island 
is — ne, with n an integer, the energy cost to add (+N, in) 
or remove ( —N , out) N extra-electrons to or from the central 
normal-metal electrode is given by 

E^\n,N) = E&'Jti ±N)~ E° wd (n) = 

= E c N 2 ±^eVN±2E c (n-n g )N, (1) 

if the tunneling process occurs through the drain (D), and 

E^\ n ,N) = E^ nd (n±N) —£f sland («) = 

= E C N 2 T \eVN ± 2 E c (n - n g )N , (2) 

when the insulating barrier of the source (S) is overcome. In 
Eqs. <[Qi and (0, the total energy of the island EF^An ± /V), 
with i = S,D, is the difference between the electrostatic energy 
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FIG. 2. (Color online) Close view of the overlapping region between the Coulomb diamonds for n = 0 and n = 1 obtained using (a) Eq/A = 1, 
(b) Ec/A = 0.6, and (c) Ec/A = 2. Also shown are the single- (dashed blue lines) and two-particle (dashed red lines) thresholds and the 
optimal loop (solid black lines) at eV ~ A from n g = n g \ to n g = n g 2 - 


due to the Coulomb interactions involving also the induced 
charge, and the work done by all the voltage sources to in¬ 
crease or decrease n with the tunneling of N particles through 
one of the insulating layers; n g = C g V g /e is the gate-induced 
charge d— 


A. Single-Electron Tunneling 

Due to the energy gap in the BCS density of states of a 
superconductor, single-electron tunneling events (N = 1) are 
energetically allowed above the gap, i.e., when the changes 
in energy Eqs. 03 and 0 are smaller than —A. On the con¬ 
trary, above —A the excess charge — ne of the island remains 
fixed to its initial value. Per each n, the threshold condi¬ 
tions £)j /out (n. 1) = —A and £( n/olU («, 1) = —A give rise to 
four crossing lines in the plot of the total bias voltage V as 
a function of the gate-induced charge n g . The four intersec¬ 
tion points between these lines are the edges of the so-called 
Coulomb diamond which is a stability region for the system. 
This means that no single-electron tunneling process can oc¬ 
cur for the values of V and n g within its area. Unlike the case 
of a fully normal SET, NININ, the Coulomb diamonds for a 
SINIS device corresponding to different n overlap. Specifi¬ 
cally, when Ec ~ A, the stability region for a given n shares 
two distinct portions of the V vs n g plane with the n + 1 and 
n — 1 diamonds, i.e., in each overlapping area at most two dif¬ 
ferent values of n are stable. This feature is at the basis of the 
generation of a controlled and synchronized single-electron 
current through the hybrid single-island structure of Fig. □ In 
this regard, let us consider, for instance, the plot of Fig. |2(a)| 
where a close view of the Coulomb diamonds corresponding 
to n = 0 and n= 1 and their shared part are shown. In princi¬ 
ple, to have a cycle corresponding to a single-particle transfer 
from the source to the drain, n g has to move along a closed 
path in the V vs n„ plane which connects the diamonds where 
n = 0 and n = 1 are stable. Thanks to the presence of the 
overlapping region, this kind of connection can be realized 
avoiding the part of the plane where both n = 0 and n = 1 
are unstable. As a result, each single-electron tunneling event 
to/from the central island can be controlled by means of the 


gate potential V g . During each cycle of n g along the working 
loop, the bias voltage V is usually kept fixed close to A/e. For 
this optimal value, the superconducting energy gap A guaran¬ 
tees an efficient suppression of thermally-activated tunneling 
events and quasi-particle excitations as well as elastic and in¬ 
elastic cotunneling processes^- 1 A typical loop used in real 
experiments with these features is shown in Fig. |2(a)| Starting 
from n g = n g i, the number of excess electrons localized on 
island, whose initial value is n = 0, remains constant until the 
threshold Si n (0, 1), defined by the equation £!/((). 1) = —A, is 
crossed. At that point one electron can enter in the central 
electrode via the source junction and n passes from 0 to 1. 
Once n g = n g 2 is reached, the closed path is covered back¬ 
ward. The extra electron on the island can tunnel out through 
the drain only after overcoming the threshold Z3 0Ut (1,1), given 
by the equation £™ t ( 1,1) = —A. When n g is again equal to 
n g i, the island is back in its initial state and a new cycle can 
start. Since per each cycle exactly one electron is transferred 
from the source to the drain, driving n g from n g ,\ to n g 2 and 
back to n g \ with a signal with frequency / allows to generate 
the single-electron current I = ef. 


B. Andreev Reflection and Higher-Order Processes 

In addition to single-electron tunneling events, the current 
flowing through a SINIS transistor is, in general, also affected 
by the Andreev reflection, i.e., the transfer of two-electrons 
per unit of time inside or outside the island . 1819 This second- 
order tunneling process is insensitive to the energy barrier pro¬ 
vided by the superconducting gap. This means that the rate of 
the transitions n —> n ±2 can be relevant although the device is 
working at the optimal bias, eV ~ A. As a result, the Coulomb 
diamonds for the Andreev reflection events are obtained just 
imposing that the energies Eqs. 0 and (0 for N = 2 are 
smaller than zero. However, as shown in Figs. |2(b j| and |2(c j] 
the energy A, together with the charging energy Ec, plays an 
important role in the determination of the two-electron tun¬ 
neling probability. If the ratio Ec/A is smaller than 1, the 
Andreev diamonds are contained within the single-particle 
stability regions. In this case, we see from Fig. |2(b)| that 
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the optimal loop crosses the two-particle threshold .S'j n (0,2), 
given by the equation £“(0,2) = 0, before the single-electron 
line 5’i n (0,1), while going from n s ,i to n g ^- When n g is de¬ 
creased back to n g _\, the closed path overcomes Z) out (l, 1) af¬ 
ter D,| Ut ( 1.2), the line corresponding to E™ 1 ( 1,2) = 0. It fol¬ 
lows that, in this regime, the control of single-electron tun¬ 
neling is compromised by the Andreev transitions 0 —> 2 and 
1 —y — 1. On the other hand, when EqI A > 1, the single¬ 
particle diamonds are smaller than the ones for Andreev re¬ 
flection. Now, the two-particle thresholds can be avoided, as 
shown in Fig. |2(c)| thereby suppressing the probability to in¬ 
crease/decrease the charge of the island by two electrons per 
each tunneling event. 

However, higher-order processes, such as the cotunnel¬ 
ing of one electron and one Cooper-pairri can occur while 
n g covers the loop of Fig. |2(c)| They can limit the single¬ 
electron transfer accuracy even if Ec/A > 1. In particular, 
the more the system stays in the overlapping region where 
more than one charge state is stable, the bigger the effect of 
unwanted transitions would be. To decrease the influence of 
the higher-order error events, the signal n g (t) which is usu¬ 
ally used to go from n g \ to n g 2 and back to n g \ is a square- 
wave. This choice guarantees that the time spent in between 
n g \ and n g 2 is minimized. On the other hand, the period 
T = 1// of n g (t ) has to be long enough in order for the single¬ 
particle tunneling processes to take place. If the number n 
changes by one electron with the rate Fi e , then the tunnel¬ 
ing error or probability that the charge of the island remains 
the same is £ ~ expf—F \ e /2f). In particular, the requirement 
e < £ me tr = 1(T 8 has to be satisfied for the definition of the 
quantum current standard. This means that T| e > 10 9 s" 1 be¬ 
cause the trade-off between the missed tunneling discussed 
above and the leakage by Cooper pair - electron cotunneling 
limits the maximum operation frequency to / ~ 60 MHz to 
achieve the metrologically-accurate current I = ef ~ 10 pA 
for a single SINIS turn stile . 1 l - 2a2i 


III. ENVIRONMENT-ASSISTED ANDREEV REFLECTION 

A. The effect of the electromagnetic environment on the 
electronic transport 

As discussed in the previous section, the tunneling pro¬ 
cesses involving more than one electron may be reduced bi¬ 
asing the SINIS turnstile at the optimal value eV ~ A, consid¬ 
ering Ec/A> 1 and using for n g (t) a square-wave-like signal 
which oscillates with frequency / between the two induced 
gate charges n g j and n g 2 of Fig. |2(c)| Under these conditions, 
one expects to measure the current I = ef with a relatively 
high accuracy. In principle, it should be possible even going 
below the relative error £ metr required by the metrological ap¬ 
plications. However, in real experiments, the achievement of 
the accuracy needed for the definition of the quantum current 
standard still remains a difficult task. 

The coupling of the hybrid turnstile with its surround¬ 
ing high-temperature electromagnetic environment may be a 
detrimental source of error— Indeed, the absorption/emission 


of energy from/to such a thermal bath allows the tunneling of 
electrons, even when the overcoming of the insulating barrier 
results to be energetically forbidden for a well isolated SET. 
Nevertheless, the environment-assisted tunneling of quasi¬ 
particles can be efficiently suppressed using, for instance, an 
on chip capacitively coupled ground planed and/or by means 
of a highly-resistive transmission line.— The main contribu¬ 
tion to the leakage current observed in the I-V characteristic is 
typically due to the Andreev reflection. Although large charg¬ 
ing energies, Eq > A, should reduce the probability for this 
two-particle process to occur, the tunneling of Cooper-pairs 
still can have a strong influence on the current flowing through 
the transistor. 12 20 The enhancement of the Andreev tunneling 
events due to the coupling of the system with the external bath 
may account for this behavior. 15 To understand under which 
conditions the environment-assisted Andreev reflection can be 
relevant, we consider the circuit of Fig. [3] where we introduce 
the effective impedances Zi(ft)), Zi (ft)), and Z„(oj) to model 
the thermal bath. We assume also that the two junctions in the 
system have the same tunnel resistance Rj. 



FIG. 3. Circuit representation of the hybrid S-I-N-I-S single-electron 
transistor (SET). The two NIS junctions constituting the device have 
the same capacitance C and tunnel resistance Rj and are connected 
to the source V$ = —V/2 and drain Vq = V/2 voltages via the 
impedances Z\ (to) and Z^iof respectively. The normal metal island 
is controlled by means of the gate voltage V g via the capacitance C g . 
The gate impedance Z g (co) together with Z\ (ft)) and Z 2 (fi>) represent 
the electromagnetic environment at temperature T ew/ . 


B. Single-photon-assisted two-electron tunneling rate 

In order to find the tunneling rate of the Andreev reflection 
process under the effect of the electromagnetic environment, 
we start by considering the tunneling Hamiltonian 

Ht = e ^ cnv ^ hep (wpip j( j T v p 7_p _(t) Qk,o + H.c., (3) 

k.p,(T 

which accounts for the transfer of two electrons between the 
normal-metal island and one of the superconducting elec¬ 
trodes of the SINIS SET of Fig. [3 Equation (0 is written in 
terms of the creation fp a (ci k ) and annihilation fp.a (^k.a) 
operators of quasiparticles (electrons) in S (N) with wave vec¬ 
tor p (k) and spin a =t,4~ The tunnel matrix elements r^.p, 
in general, depend on p and k. The BCS coherence factors u p 
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and v p are spin-independent and satisfy the relations 


2 i 2 

w p 1 Pp 


1 

2 



A 

" PVP _ ( £ 2-A2)l/2 


(4) 


where <i; p is the energy of an electron in S with momen¬ 
tum p measured with respect to the Fermi level, and e p = 
(c 2 + A 2 ) 1 / 2 is the quasiparticle energy. The translation op¬ 
erator e'fe™ in Eq. 0 accounts for the change of the charge 
of the electrodes due to the environment-assisted tunneling of 
one electron. Considering the environment as an infinite en¬ 
semble of quantum harmonic oscillators with temperature T env 
(Caldeira-Leggett model22 - ™), the fluctuating phase <p env can 
be written as 

<Penv =TjP^ (4+^0 (5) 

2. A 


where the phase <p x represents the position operator of 
the harmonic oscillator A with mass C- A and characteris¬ 
tic frequency 0 J;_ = 1 /y/LxCx- The coupling term is px = 
(e/Ti)y/h/2Ci(Oh, and the operators c/ and £x create and an¬ 
nihilate one photon with energy JUx>x (see Appendix^- Here¬ 
after, we assume that the coupling of the SINIS with the envi¬ 
ronment is weak, meaning that at most a single photon is in¬ 
volved in the exchange of energy between the system and the 
thermal bath: 14 - In other words, we consider the limit where 
Px •Cl and the series expansion of the charge translation op¬ 
erator Eq. (0) in Hj can be truncated at the first order, i.e., 
e i<Pi:m ~ i _|_ /(p env . The validity of this assumption will be dis¬ 
cussed in the following. 

Let us focus on the Andreev process 1 —> — 1, character¬ 
ized by the transfer of two electrons from the normal metal 
island to the superconducting drain electrode as a Cooper pair. 
According to perturbation theory in //y , the total probability 
amplitude to have such a second-order event in the system of 
Fig0is given by 



FIG. 4. Feynman diagrams of the two possible single-photon ab¬ 
sorption processes giving rise to the environment-assisted Andreev 
reflection. 


one of the two tunneling electrons can absorb the energy of the 
only available photon, in the weak coupling limit. As a result, 
for a fixed wave vector p of the virtual quasiparticle in S, only 
the four intermediate states 


11 a) — l k l T® Impairs , p)s ® «A 4“ l)env , 

I^a) = l k l t® Impairs, p)s ® | H A)env > 

|3a> = [£*<£, k 2 4 -)#® I ^pairs 7 p) S ® K + i) env i 
Ka) = U&Kt/k? A)*'® |UpairsiP)s® | n A)env j (9) 

can give a non-zero contribution to A£ k . The difference 
between the energies of these virtual states and the en¬ 
ergy — if] of the initial state \ix) determine the ampli¬ 
tude Eq. (0). The imaginary part r) = hY\ >o/2 accounts for 
the lifetime broadening of \ix) due to the competing single¬ 
electron tunneling processes occurring with rate r^o- Ac¬ 
cording to perturbation theory in the tunneling Hamiltonian 
//y , the first-order rate, describing one electron going out of 
the island through the drain, can be written as 


yDynes _ 

1 i->o — 


1 A R k f\ E D l YA)\N^ nc \E/ A) 


J_A Rk f 

2 K h Rj Jo 


dE (10) 


in terms of the Dynes density of states of a superconductor, 25 


A {fx\HT\m x ) (mx\H T \ix) 

k i.k 2 L r r. _i_ ; r 


(6) 


N^(E/A) 


91e 


E/A + fy Dynes 

VWA + IXJ 2 -! 


( 11 ) 


for fixed values of the environmental index A, and of the initial 
wave vectors ki and k 2 . Here the initial state is 

la) = l k l t,k 2 A)iV® |«pairs, |«A + l)env , (7) 

with two electrons in N with opposite spin and momenta ki 
and k 2 , « pa irs Cooper pairs and no quasiparticle excitations in 
S, and nx + 1 photons with energy ha>x in the environment. 
On the other hand, the final state is 

I/a) = <8 Impairs + 1 , <8 |«;i)env , (8) 

with an additional Cooper pair in S, two less electrons in N, 
and one less photon in the Caldeira-Leggett bath. The tran¬ 
sition from the state (0 to the state 0 is determined by all 
the possible intermediate virtual states | mx) such that a quasi¬ 
particle with momentum p is created in S after the annihilation 
of one of the two electrons in N. As illustrated in Fig. 0, only 


which depends on the phenomenological Dynes parameter 
y Dynes . In Eq. ®, E° ut (l, 1) = 2 E c (n g - 1/2) - eV/2 is the 
energy cost that has to be payed by the voltage sources in or¬ 
der for the transition 1 —► 0 to occur [see Eq. 0| -R K =h/e 2 
is the resistance quantum. The Dynes rate Eq. (ITOt is valid in 
the zero-temperature limit, AbTsinis "C A, and takes into ac¬ 
count the most relevant single-electron error sources, such as 
the environment. 

Using Eqs. 0. (0, and 0. the amplitude Eq. 0 reads as 
A k!,k 2 = i to VfkiVM PAV^I E ( m p v p) 5 pA > 0-2) 

p 

for a low-temperature hybrid single-electron transistor, 
ksT sinis A, and assuming constant tunneling matrix ele¬ 
ments, fk.p = ?k p = to (point tunnel junction). In Eq. (fl2l i. 
we introduced the Fermi-Dirac distribution function /k = 
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FIG. 5. (Color online) Photon-assisted Andreev rates, given by the numerical evaluation of Eq. oi. as a function of the gate-induced charge 
n s with A = 210/leV (Aluminum), Rj = 430 k£2, C s = 10~ 16 , JV = 100, and )fo yncs = 10~ 5 . In panel (a), for each rate R = 1100£2 and 
Ec/A =1.4 with C = 0.86 • 10 16 ; the values of r env are: 70 mK (red), 140 mK (blue), 780 mK (green), 1.5 K (orange), 4.2 K (purple). In 
panel (b), for each curve T env = 1.5 K and Ec/A = 1.4 with C = 0.86- 10~ 16 ; the resistances R are: 1100 £2 (red), 10 £2 (blue), 0.1 £2 (green), 
0.001 £2 (orange). In both panel (a) and (b), the dashed black line is the Andreev rate valid in the absence of environment (see Ref. 11). In 
panel (c), for fixed R = 1100 £2, the curves with the same color are obtained using the same charging energy, Ec/A: 1.4 with C = 0.86 ■ 10~ 16 
(red lines), 1.8 with C = 0.558 • 10~ 16 (blue lines), 2.5 with C = 0.262 ■ 10~ 16 (green lines); the values of T env are: 4.2 K (solid curves), 
500 mK (dashed curves), and 100 mK (dotted-dashed curves). In all the three panels, also shown are the single- and two-particle thresholds, 
1/2 — A/4 Ec (light-blue vertical dotted lines), and A/4£/■ (light-red vertical dotted lines), respectively. 


[exp((i;k/£g7siNis) + 1] 1 for the normal metal electrons and 
the sum of the intermediate-state denominators 


Sp,A = 


+ 


1 

£p-<?k, + irj 
1 

£p - <?k 2 + ir\ 


1 

e p - <?k 2 - ha)x + ill 
1 

-h(o x +iri 


(13) 


Here £ p = e p + E/f( 1,1) is the virtual state energy and <g k is 
the energy of an electron in N with momentum k measured 
with respect to the Fermi level. Summing over all the possible 
initial states and considering the spin degeneracy, one obtains 
the total rate 


penv _ V V" 1 

1 ar - -r L L 


kj.ki A 


1 ki ,k2 


s(^ lM2 +n(o x ) , (14) 


where kz = Ck| + q,ki — Ep" (1,2) is determined by the en¬ 
ergy cost E// l ( 1.2) —4 Ecn g — eV needed for the second-order 
transition 1—2—1 to occur [see Eq. ©]■ The environment- 
assisted Andreev rate Eq. (IT4l) is written in terms of the prob¬ 
ability 


a2 

A k!,k 2 


2 


to fkjk 2 Pi n x £ ( m p v p) (w) 5 p ,a 5 p',a • 
p,p' 


Approximating the sums over ki, k 2 , p and p' by the corre¬ 
sponding integrals, assuming that n x is given by the Bose- 
Einstein distribution, hbe{(0x) = [exp/Tzco^/ kgT env ) — 1] , 

and using the properties of the Dirac delta function, Eq. ([T4l i 


can be written as 


-p env 
1 AR — 


X 


(MpVp) (wp'Vp/) fkjk 2 >ibe{-£, k 1 ,k 2 )/(-4k 1 ,k 2 ) 


x ( 5 P 5 P') LPa^S^^+^Oa) • (15) 

A 


Here S p and S*, are .S' p ;_ and S*, ^ evaluated for ho) x = 
—| k k . Thanks to fluctuation-dissipation theorem, the sum 
over A in Eq. (IT5l) can be related to the effective impedance 
9ic [Zeff(ffl)] “seen” by the drain (see Appendices [A] [B] and 
©. Then, in the low-temperature limit A:b7sinis < A, Eq. (d 
can be recast in the form 

A 2 n 2 


t -1 env 
1 AR — 


1 


{/Inf TiR k JV 


r k \ f° f + °° 


x (V £ P- A2 \Ap'” a2 ) nBE (-^,k 2 ) / Hk,.kJ 

x (s p s^) 91e[Z eff (4 c ljk2 //l)] , (16) 


using the BCS relation for w p v p given in Eq. ©, and the 
quasiparticle energies £ p and £ p / as integration variables. In 
this last formula, we also introduced the number of conduct¬ 
ing channels of the junction. The environment-assisted 
Andreev rate Eq. (IT6t is valid in the single-photon regime 
px <C 1, i.e., for small values of the effective impedance, 
91e [Z e ff(ffl)] /Rk -C A/kBT enY .— Additionally, since we as¬ 
sumed that only the absorption process can occur, one has to 
impose that | k k < Pin Eq. dTTb. namely Ej( ut ( 1.2) > 0. This 
means that Eq. ( fl6l i applies only for those values of the in¬ 
duced gate charge n ? equal to and larger than the two-particle 
threshold eV/4Ec- 
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(a) (b) 




(c) (d) 

FIG. 6. (Color online) Plot of the ratio e acc as a function of the gate-induced charge n g with A = 210 /leV (Aluminum), Rj = 430 k£2, 
C g = 10~ 16 , and JV = 100. In panel (a), for each rate y Dy „ s = 10 5 , R = 10 £2 and Ec/A =1.4 with C = 0.86 ■ 10~ 16 ; the values of r e nv 
are: 70 mK (red), 140 mK (blue), 780 mK (green), 1.5 K (orange), 4.2 K (purple). In panel (b), for each curve yy yn( , s = 10~ 5 , T env = 1.5 K and 
EqI A= 1.4 with C = 0.86 • 10~ 16 ; the resistances R are: 110012 (red), 10 £2 (blue), 0.1 £2 (green), 0.001 12 (orange). In panel (c), for fixed 
Tbynes = 10~ 5 and R = 10 £2, the curves with the same color are obtained using the same charging energy, Ec/A: 1.4 with C = 0.86- 10~ 16 
(red lines), 1.8 with C = 0.558 ■ 10~ 16 (blue lines), 2.5 with C = 0.262 • 10~ 16 (green lines); the values of T env are: 1.5 K (solid curves), and 
500 mK (dashed curves). In panel (d), the Dynes parameter y Dyncs is equal to 10~ 4 (solid lines) and 10~ 7 (dashed lines). The curves with the 
same color are obtained using the same charging energy, Ec/A: 1.4 (red lines), 1.8 (blue lines), 2.5 (green lines). In all the four panels, the 
single- and two-particle thresholds, 1/2 — A/AEc (light blue vertical dotted lines), and A/4 Ec (light red vertical dotted lines) respectively are 
also shown. 


For the circuit depicted in Fig. [3j the [Z e ff(a>)] is equal to 
the right-hand side of Eq. (IB4b . The latter is the sum of three 
terms which are of the same order of magnitude for the typical 
experimental values of the capacitances C and C„ ~ C, and 
of the impedances Z\{(o) ~ Zn{(o) ~ Z g (a>). Consequently, 
because of the symmetry of the circuit of Fig. [3] with respect 
to the gate, we focus hereafter on the case where the voltage 
fluctuation across Co is produced only by Z g (co), considering 
Zi(ffl) and Z 2 (ft>) as noiseless. In addition, we assume that 
the effect of Z\ and Z 2 can be neglected with respect to Cs 
and Cd, namely Z 12 [£//*(I -2)/ft] <C h/E™ 1 (\,2)Cs t D- As a 
result, setting (5V"i = 5V 2 = 0 as well as Z\ (w) = Z 2 (co) = 0, 
Eq. (1B4I) yields 


91e [Zeff(ft))] 


R 

(C z /C s y- + (2coRC) 2 ’ 


(17) 


for a purely resistive environment Z g (co) =R and a symmetric 
turnstile Zq. = Zq d = Zq. 


In the particular case where 91e [Z e g-(o)] can be approxi¬ 
mated with a frequency independent resistance R , Eq. ( 1 1 61 ) 
becomes 


-hT , 
AR ' 


/env 


(2 n) A 




x ' (* i .^)/® lfc ) 2 . as ) 


in the high-temperature limit, ksT env >^ ut (l,2)=4£ c ^- 
A with (1 /4) < n g < (3/4), assuming that the system is work¬ 
ing at the optimal point eV = A, and for large charging en¬ 
ergy Ec > A. In Eq. ( ITSl ). we introduced the high-temperature 
Dynes parameter = 2k(R/ RK){kBT em /A) which is 
the only term of f which depends on R and r env . 
































7 


C. Numerical Results 

Using Eq. (ITTl) . the numerical integration of Eq. ( 1 1 6[ i is rel¬ 
atively straightforward. Figure [5(a)| shows the photon-assisted 
Andreev rate, Eq. ( fl6l ). as a function of the gate-induced 
charge n g , for a single-electron transistor biased at the opti¬ 
mal voltage, eV = A, and with charging energy Eq > A. Each 
curve is obtained for different values of the temperature of the 
environment T env . The other parameters are fixed to the val¬ 
ues of sample S3 of Ref. [u|, as indicated in the figure. We 
see that the probability to have the tunneling of a Cooper-pair 
can be different from zero also away from the two-particle 
tunneling threshold, unlike the case without environment. In 
particular, the exchange of energy with the thermal bath in 
which the SET is embedded can make the Andreev reflection 
relevant even around the single-particle threshold. As a result, 
although the boundary of the Coulomb diamond correspond¬ 
ing to the transition 1 —> — 1 is avoided by means of the loop 
of Fig. |2(c)| a Cooper pair can tunnel through the barrier of the 
drain, while n g goes back to n g \, before crossing the 1 —> 0 
line. The decrease of T env leads to smaller values of [see 
Fig- [5(a)), as well as the use of an electromagnetic environ¬ 
ment with a smaller resistance R [see Fig. |5(b)[ . Whereas, in 
the latter case, the whole Andreev rate curve is shifted down 
proportionally to the ratio between initial and final resistances, 
the modulus of the first derivative of Eq. ( IT6t for n g > A/4 Eq 
increases proportionally to Tenv [see Fig. |5(a)[ . 

The dependence of the photon-assisted Andreev rate, 
Eq. (1 1 6b . on the charging energy Ec is shown in Fig. |5(c)| 
The increase of the ratio Ec /A > 1 allows to reduce the effect 
of the two-particle tunneling on the total electric current sus¬ 
tained by the SINIS turnstile. In particular, the lower is the 
effective temperature of the environment with respect to the 
critical temperature of the superconductor, the larger is the re¬ 
duction of r™^ upon increasing Eq/ A. Notice that the main 
effect of the change of the charging energy Eq is the shifting 
of the environment-assisted Andreev rate along the induced- 
gate charge axis by the difference between the initial and final 
inverse ratios A/4 £/-. 

Assuming that the number of electrons of the metallic is¬ 
land of the circuit of Fig.[3]decreases because of the tunneling 
of quasi-particles and Cooper pairs only, the total rate can be 
written as 

T~' T-'Dynes . 

1 tot — i l—^o + AR ■ 

As a result, the error e aC c = ^F^/Fj/^q determines how much 
the environment-assisted Andreev reflection affects the charge 
transport in the SINIS transistor. In particular, the condi¬ 
tion e acc < 10 x is required for the metrological applications.— 
Figure ([6ji shows the ratio £ acc obtained from a numerical 
evaluation of Eqs. (fl~6t and ( ITOt . as a function of n g , when 
Eq. ® holds. We see that £ acc is a non-monotonic func¬ 
tion of rig. Starting from the two-particle threshold occur¬ 
ring for rig = A/4E C , this error first decreases exponentially 
as n g is increased. Then, close to the single-particle thresh¬ 
old, it rises up again reaching a local maximum value around 
n g = 1/2 — A/4 Ec- For larger n g it tends exponentially to 
zero. Because of this kind of behavior, £ acc can be smaller 


or of the order of 10 -8 when A/4 Ec < n g < 1/2 — A/4 Ec, 
and, at the same time, much larger than the value required 
by metrology around the single-particle threshold. Conse¬ 
quently, the time spent by the signal used to drive n g around 
n g = 1/2 — A/4 Ec has to be as small as possible in order to 
minimize the environment-assisted Andreev reflection. 

From the experimental point of view, the determination, 
with a relatively high accuracy, of the values of the effective 
parameters of the environment, R and Tenv, is a tough chal¬ 
lenge. The use of the Dynes parameter y Dynes , which in general 
depends also on R and T env , is preferred because it can be de¬ 
termined from the measured current-voltage characteristic of 
the SINIS turnstile. In this regard, the high-temperature two- 
particle tunneling rate Eq. ( IT8l > allows to study the photon- 
assisted Andreev reflection in terms of y Dynt , s only. In Fig. |6(d)| 
we plot the error £ acc obtained using Eq. (ITSl) as a function 
of n g . We see that the Dynes parameter, which typically 
ranges from 10 4 to 10 7 , strongly affects r/j 1 , in the range 
A/4 Ec <n g < 1/2 —A/4 Ec- On the contrary, y Dynes plays a 
minor role in the reduction of £ acc when n g is close to the 
single-particle threshold. 



FIG. 7. Circuit representation of the hybrid S-I-N-I-S single-electron 
transistor (SET) connected to the impedances of the electromagnetic 
environment Zi(co), Z 2 (co), and Z g (co) by means of three transmis¬ 
sion lines, (1), (2), and (g), respectively. 


IV. EFFECT OF A RESISTIVE TRANSMISSION LINE ON 
THE PHOTON-ASSISTED ANDREEV RATE 

The results presented in the previous section have been ob¬ 
tained considering a SINIS turnstile directly connected to the 
external electromagnetic environment as illustrated in Fig. [3] 
Using Eq. O, we have shown that the smaller are the tem¬ 
perature 7en V and the resistance R of the external circuit, as 
well as the total capacitance Cz, the lower is the Andreev tun¬ 
neling rate Eq. (IT6t with respect to the single-particle one. In 
particular, we have seen that the metrological accuracy may 
be reached for certain values of r env , R, and Cz- However, in 
real experiments, the control of these parameters is typically 
limited. In general, their tuning to the desired values can be 
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FIG. 8. (Color online) (a) Plot of the photon-assisted Andreev rate Eq. {Hi and (b) of the corresponding ratio e acc as a function of the gate- 
induced charge n g . In panel (a), the solid and the dashed lines are obtained using Eqs. {H and CEO. respectively, setting zo = c o = r o = 0 
(red lines) and using zo = 0.7 and Co = 1 with ro = 5 ■ 10 3 (blue lines), ro = 5 ■ 10 4 (green lines), and ro = 5 • 10 5 (orange lines). In panel (b), 
the parameters for the plotted ratios are zo = co = ro = 0 (solid lines) and zo = 0.7 and co = 1 with ro = 5 ■ 10 4 (dashed lines) and ro = 5 ■ 10 5 
(dotted-dashed lines). The curves with the same color are obtained using the same charging energy, Eq/A : 1.4 with C = 0.86- 10~ 16 (red 
lines), 1.8 with C = 0.558 • 10~ 16 (blue lines), 2.5 with C = 0.262 ■ 10~ 16 (green lines). In both panels, we set A = 210 peV (Aluminum), 
Rt = 430 k£2, C g = 10 16 , ,jY = 100, y Dyncs = 1(G 5 , T e nv = 1.5 K, and R = 10 £2. The single- and two-particle thresholds, 1/2 —A/4 Eq 
(light-blue vertical dotted lines), and A/4£/- (light-red vertical dotted lines), respectively, are also shown. 


a difficult task. As discussed in Ref. 0 the insertion of cold 
and lossy transmission lines between the turnstile and the en¬ 
vironment can help in overcoming this problem. One expects 
that such an indirect coupling allows a further reduction of the 
environment-assisted two-particle tunneling. 


We therefore consider the circuit of Fig. [7] where the three 
impedances of the environment are connected to the SINIS 
turnstile by means of three transmission lines. We assume 
that the latter are noiseless, i.e., at zero temperature. The 
noise across the drain of a SINIS device in such an indirect 
configuration is derived in Annendixlcl In this case, the effec¬ 
tive impedance 9Ie |Z e n/oj)] appearing in Eq. ( 1 1 61 ) is given by 
Eq. (lC6b . Let us assume again that Z\(co) ~ Zn{(o) ~ Z g (co) 
and C ~ C g as well as that the three transmission lines have 
length l and are all described by the same parameters Ro, 
Co, and Lq, the resistance, the capacitance, and the induc¬ 
tance per unit length, respectively. Then, the three terms in 
the right-hand side of Eq. (lC6b contribute in a similar way to 
9Ie [Zeff(fi))]. In particular, they have the same order of mag¬ 
nitude for large Rq and l. On the basis of these considerations 
and given the symmetry of the circuit of Fig. [7J we assume 
that the detrimental noise comes only from Z g {(o), i.e., the 
voltage noises are 8V S ^ 0 and 8V\ = SVi = 0, neglecting the 
effect of Z\ (ft)) and Z 2 (co) as well as of the transmission lines 
(1) and (2). As a result, the effective impedance “seen” by the 
drain reduces to 


9Ie[Z eff (©)]«R M®)[ 2 \T S ((0)\ 2 , (19) 


with 


a g (co) = - 


+ ^ (Ac * + 1) 


+ An 


Z g +1 


0 i K s(m)t 


-1 


T g (co) 


( 20 ) 


for Z g (co) = R and under the condition Z\ 2\E^ l (\,2)/K\ -C 
Tl/E™ 1 {\,2)Cs,d- The effective impedance Eq. ( | 1 9b tends to 
the asymptotic expression 


9Ie [Zeff(co)] ~ R 


C \ 2 e -V 2(Bfi oCo 
2 CJ 1 + coR 0 C g /Co ’ 


( 21 ) 


if the transmission line is highly resistive, Ro 3> 
LoE™ x (\,2)/Ti, and long enough, £^/2E^ lt (l,2)RoCo/h 3 > 1 , 
and when the resistance of the environment is small, 
R <C R()Cg/2Q)J-^ We see that Eq. (f2TT) and in turn the 
environment-assisted Andreev rate Eq. ( IT6b decay exponen¬ 
tially in terms of £ and Ra¬ 
in Fig. [8j we show the plots of the Andreev rates and the 
ratios e acc resulting from the numerical integration of Eq. ( 1 1 6| | 
with 9Ie IZeiy(oj)] given by Eqs. ( IT9t and (f2Tl) . In both panels, 
we used the dimensionless parameters zo = y/To/Co/R, cq = 
£Co/C, and ro = £Rq/R whose values are chosen according to 
the analysis about the transmission function given in Ref. 0 
and in agreement with the currently achievable experimental 
values of £, Lq, Cq, and Rq. Similar results hold also in the 
configuration obtained restoring the fluctuations due to Z\ (ft)) 
and Z 2 {co) and treating Z,, (ft)) as a noiseless impedance. 

From Fig. |8(a)| we see that the bigger is ro the smaller is 
r^. In other words, a long and highly resistive transmission 
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line allows a reduction of the environment-assisted Cooper- 
pair tunneling. By comparing the solid and the dashed lines 
obtained using Eqs. m and <Efl- respectively, we note that 
this decreasing is exponential-like. Eventually, the increase of 
i and/or Rq leads to a decrease of e acc below 10 8 , even close 
to the single-particle threshold [see Fig. |8(b)| . As a result, the 
use of a highly resistive and noiseless transmission line allows 
to filter out effectively the photon-assisted Andreev tunneling 
and, in particular, to reach the accuracy needed for metrologi¬ 
cal applications. 


V. CONCLUSIONS 

In this paper, we studied the environment-assisted Cooper- 
pair tunneling in a SINIS turnstile working in the Coulomb 
blockade regime. Specifically, we derived the Andreev re¬ 
flection rate when only a single photon of the thermal bath 
is involved in the process. We found that the single-photon 
absorption enhances the two-electron tunneling from N to S. 
In particular, the probability per unit of time to have An¬ 
dreev events is different from zero even for values of the 
induced gate charge n g close to the single-particle threshold 
1/2 — A/4 Ec- As a result, the single-electron current, which 
is expected to be the dominant one in the device when n g fol¬ 
lows the loop shown in Fig. |2(c)[ is also affected by the tun¬ 
neling of Cooper pairs due to the environment. The influence 
of this source of error on the total current can be reduced by 
decreasing the effective resistance R and temperature r env of 
the environment or, equivalently, the Dynes parameter y Dyncs . 
The achievement of the metrological accuracy is also possible 
with the increasing of the charging energy Ec with respect to 
the superconducting energy gap A. We finally show that using 
a cold and lossy transmission line to couple indirectly the en¬ 
vironment with the SINIS turnstile allows to reduce further the 
probability to have two-electron tunneling events, especially 
when n g crosses the single-particle threshold while covering 
the optimal loop. 
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Appendix A: Caldeira-Leggett model and 
Fluctuation-Dissipation Theorem 

According to the Caldeira-Leggett model, the impedance 
Z(ft)) of an electric circuit can be modeled as an ensemble of 


infinite quantum harmonic LC oscillators with Hamiltonian 


He 


E 


@X I 1 r „2 

2Q + 2 Q ^ 



The charge Q- A and phase (hj e)(p- L operators play the role 
of the momentum and position respectively of the parti¬ 
cle/oscillator A with mass C; and characteristic frequency 
coj = 1/L A C A . Each oscillator A of the ensemble/environment 
affects both the charge Q and phase (p of the circuit. In partic¬ 
ular, the total phase fluctuation <p env of (p due to Z('ft)j is given 
by the superposition of all the phases of the oscillators of the 
environment, i.e., <p env = Y.X<Px- Since <p A is the position op¬ 
erator of a harmonic oscillator, <p env can be written as 

*Penv = Y. Pa ( 'c\ > (Al) 

X 


in terms of the creation c| and annihilation c A operators of 
one photon. In Eq. ( 1A1I) . we introduced the coupling term 
Px = (e//j)y / /i/2Qo)^. In the Heisenberg picture, <p env de¬ 
pends explicitly on time, with ct (t) = e +, ®V ct and c A (/) = 
e im> - t c~ A . The first time-derivative of Eq. (I A1 b gives the fluc¬ 
tuating voltage operator 


Venv(0 


h d <Penv (0 

e dt 


Ti 

e 


EPA 

X 



exit) 


(A2) 


whose mean value over the eigenstates of H env is zero. 
On the other hand, the voltage-voltage correlation function 

SVenvM) = <{Venv(/),Venv(0)}) is 

<5Venv(/,0) = f-') Y^PxPx' {in a ) i ia >X’) Q,(f,0) , (A3) 
W XX' 


with 


C A (f,0) = ({[4(/)-c A (/)],[4(0)-c A (0)]}} . (A4) 


The symbols {,} and (...) in Eq. ( lA4b indicate the anti¬ 
commutator and quantum mean value over the eigenstates of 
Henv, respectively. Assuming that the number of photons of 
the environment is infinite, the terms in Eq. (IA3b which create 
or destroy more than one photon can be neglected. Conse¬ 
quently, the correlation function (W env (f,0) becomes 


SVenv(CO) ~ 0) ZplaUSW + e-^il+tox) , 

(A5) 

where n~ A is the mean value of photons with frequency ft); . 
The Fourier transform of Eq. (1A5I) gives the spectral density 
function of the thermal bath, 

x 2n [8(co - a>x) + S(co + a>x)] ■ (A6) 
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To obtain Eq. (IA6b we assumed that nx is given by the Bose- 
Einstein distribution function = [exp(/zft)^/A:B7’env)— 

1]~* which satisfies the relation 1+2 hbe(x) =coth(x/2). T env 
is the temperature of the environment. 

On the other hand, assuming that the Fourier-transformed 
correlation function [5Venv(T0)] satisfies the quantum 
fluctuation-dissipation relation 

[5Venv(f,0)] =2ha> 9te[Z(e>)] cothfi-^A . (A7) 

\Z /cgi e nv / 

Comparing Eq. dA6b with Eq. dA7b , we finally get the expres¬ 
sion 

9*e[Z(®)] = ^-EPa®a [5(o-ft)A) + 5((0 + ft)A)] , (A8) 
z X 


hold. They can be rewritten as 

2c, SVi = (Z c , + Z,) AT, - Z,A+ c 
Z C2 SV 2 = - (Z C2 + Z 2 ) A+ 2 + Z 2 AV c 
Zc g SV g = - (Z Cg + Z g ) AVg + ZgAV c 

0 = -Z Cg Z c2 AVi-Z Cg Z cl AV 2 -Z Cl Z C2 AV g 
+ ( z c g Zc 2 +Zc g Zc l +Z Cl Zc 2 ) AV C , 

in terms of the voltage noises 8V\ = Z\8l\, SV 2 = Z 2 Sl 2 , and 
8V g = ZgSlg. By solving this system of equations one can get 
the unknown potentials A+i, AV 2 , AV S , and A+ c . After some 
algebra, the voltage drop 8V reads as 

s v = ^^i2r l {co)S v l + 2 2 {co)8v 2 -2r 3 (co)SVg \, (Bi) 


which allows to relate the macroscopic impedance Z(ffl) with 
the microscopic quantities characterizing the environment. 


Appendix B: Voltage fluctuations across the drain 

In this Appendix, we consider the circuit of Fig. [9] We 
proceed in the evaluation of the voltage noise 8V = AV C — 
AV 2 across the capacitor C 2 and of its correlation function. 
The latter can refer, for instance, to the drain of the SINIS 
transistor. 



Z 2 


where we introduced the impedances 


S(co) = ir 3 (®)[Z 2 (©)+Z C 2 (©)]/Z C2 (©) 

+ ^l(®)[JS(®)+ Z 2(®)+Z C2 (©)]/Z C2 (fi)) , 

S l (co) = Z g (co)+Z Cg (co) , 

S 2 (co) = Z\ (fit) +Zc,(fit) +Z g (fi>) +Zc g (fit) , 
^ 3 (fit)=Zi(fit)+Z Cl (fit) . (B2) 

Here Zc ; (fit) = i/(coCj ) is the impedance of the / 1 h capaci¬ 
tor with j = 1, 2.g. If 5Vi, 5+2, and 8V g satisfy the quan¬ 
tum fluctuation-dissipation theorem [see Eq. (1A7M . then the 
voltage-voltage correlation function of 5+ is 


[SV(t,0)\ a) = 2h(o 
“h 2 flCO 

2hco 


S(co) 


S’(co) 
S 2 (co) 2 


ar{(o) 

s 3 (<o) 


91e [Zi (fit)] coth 
91e[Z 2 (fit)]coth 


S’(a) 


1 hco \ 

2k^\) 


1 hco 


^[Z g ( ®)] coth - 


2 k B T 2 
1 hco 


2 kBT„ 


(B3) 


FIG. 9. Circuital scheme of a SINIS turnstile connected to an elec¬ 
tromagnetic environment which produce current noise. 


We start by considering clockwise currents in the two 
meshes of the circuit of Fig. [9] Then, assuming that the 
impedances Zi, Z 2 , and Z g produce the current noises 5/i, 
5/ 2 , and 8l g , respectively, the following equations 


A h = 8l\ - 
AI 2 = 8l 2 + 
A L = 8lo + 


A+j 

zT 

A+2 

Z 2 

A+„ 


AI g = Ah-AI 2 , 


A+c — A+i 

Zc, 

A+ c -A+ 2 
Z C2 ’ 
A+ c - A+g 


In Eq. (lB3b . we introduced the temperatures T\, 7+ and T g 
of the impedances Zi, Z 2 , and Z„, respectively. Similarly to 
Appendix lAl the function [5+(f,0)] can be related to the 
microscopic properties of the environment acting on the ca¬ 
pacitance C 2 by comparing Eqs. (lAfil ) and (IB3b : 


Rk 


EpI®a 5( |e>|-e>ji) = 


S(co) 


S((o) 

S 2 (co) 


+ 


S(<o) 
S3 (a) 


S(co) 


^e[Zi(®)] 


91e[Z 2 (®)] 


9le[Z g (fi))] . (B4) 


To obtain Eq. (1B4I) we imposed that the temperature of the 
environment is uniform, T\ = T 2 = T g = T env . 
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Appendix C: Voltage fluctuations across the drain in the 
presence of transmission lines 


In this Appendix, we analyze the circuit of Fig. 1101 where 
the transmission lines (1), (2), and (g) are inserted between 
the capacitances of the turnstile C i, C 2 , and C g and the ef¬ 
fective impedances Z\, Z 2 , and Z g which give rise to the cur¬ 
rent noises 8l \, 8h, and 8l g , respectively. The /'th line, with 
j = 1,2 ,g, has length and is characterized by the param¬ 
eters Rq \ Cq \ and Lq \ the resistance, the capacitance, and 
the inductance per unit length, respectively. As in Appendix 


Vi(0) Cl A F c C2 V 2 (0) 



0 


h 

X 


Z 2 


icoL^)/Kj((o) [see also Refsfl4landH7ll. The unknown coef¬ 
ficients Aj and Bj and the potential drop AV C are the solutions 
of the system of equations presented in Eq. (1C Ik After some 
algebra, one can derive the voltage noise SV across C 2 : 


8 V = a\(co)T\(co) 8 V\ + a g (a>)T g (co) 8 V g + ch((o)T 2 (co) 8 V 2 . 

(C3) 

Here dVi = Z\ 8 l\, 8 V 2 = Z 28 /?, and 8 V g = Z g 8 l g . In 
Eq. ( IC3b . we introduced the coefficients 


( \ 4(e>) ( ^ 1 

cq ,,(co) =-;— 0,2(60) = 1 — 

^ ' Z Cl Jco)Y((o) 2V ; 


with 


ZcAco)Y(co) ’ 
(C4) 


F(co) = 1(1 - Ac 2 ) + A 2 e iK ^T 2 (co) , 


Y(co) = lc ‘ 


1 

2zi 1) 

Aq+_1 

4) 


- X\ 


Zc { +1 
X\ +1 
+ 1 
X g + 1 


iKi(co)ti 4(<H) 

Z c ,(co) 

JK g (a>y g 4M 
Z Cg (a>) 


2 Z [ J> \z g + ij 

Zc 2 + 1 -\~X 2 ( Ac ' 2 A 1 1 e iK 2 ((0)i 2 T 2 (co) 


where 


FIG. 10. Circuital scheme of a SINIS turnstile connected to a noisy 
electromagnetic environment by means of transmission lines. 


m we focus on the derivation of the voltage-voltage correla¬ 
tion function of the potential SV = AV C — V 2 (0) across C 2 . 

Imposing clockwise currents in the two main meshes of the 
circuit of Fig. [TOj one can write the system of equations 


El (It) 

Zi 

Zg 

Vlih) 
Z 2 ’ 

4(o) = / 1 (o)-/ 2 (o), 


h(i x ) = 8h- 

4(4) = ^4 + 

h(f- 2 ) = Sl 2 - 


where 


4 ( 0 ) = - 


AV c -Vi(0) 


> ‘g 


4(0) = 


4(0) = 


Zci 
AV4-V g (0) 
Zc : 

Ay c -v 2 (o) 

Z(_2 


Vj(x) =Aje iK ^ x + Bje- iK J^ x , 


Ij(x) = 


Aje iKi ^°) x — B je ~ iKj ^ x 


(Cl) 


/Z { J\co) (C2) 


are the voltage and the current at a given point x along the /th 
transmission line. In Eqs. ( IC2k we introduced the wave vec¬ 
tor Kj(co) = 60 2 Lq' > Cq' 1 +i(OR { 0 j) ctf ) of the signal propagat¬ 
ing along the jth line, and the impedance zii\(o) = i(/?| A — 


Tx((0) = 
t 2 (co) = 
Tg(co) = 


1_ (1 — Ac t )(Ai + 1) 

2 AiA q e~ iK x ( ffl )4 — e iKl (®)^i 

£ (1 -Zc 2 )(Z 2 + 1) 

2 e -iK 2 (a>)e 2 _ A 2 A C2 e iX 2( £8 )4 

1_ (1 — kc g )(k g +1) 

2 e ~— A„Ar 

6 '-g 


are the transmission functions, and 


A j(C0) 


Z { J\(o)-Zj(a>) , 

Z { J\(o)+Zj((o) ’ 


A q (®) 


Z { J\co)-Z Cj (co) 

Z ( 2\(o) + Zcj((o) 


the reflection coefficients^ 4 

Assuming that 5Vi, SV 2 , and 8V g satisfy the fluctuation- 
dissipation theorem [see Eq. (IA7H . and thatZi, Z 2 , and Z„ are 
at the same temperature, we can finally get the relation 


^-Y J Pl^(\co\-co x ) = \ax(co)\ 2 14(a))! 2 Heft (a)] 
z x 

+ |a 2 (®)| 2 |4(®)| 2 SRe[Zj(fi))] 

+ |a 4 ,(a))| 2 |4(a))| 2 9te[Z g (fi>)] . 

(C6) 


by means of Eq. ( 1 A 6b . 
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